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ABSTRACT 

■ The theory of cold, relativistic, magnetohydrodynamic outflows is generalized by the inclusion of an intense 

radiation source. In some contexts, such the breakout of a gamma-ray burst jet from a star, the outflow is 
heated to a high temperature at a large optical depth. Eventually it becomes transparent and is pushed to a 
higher Lorentz factor by a combination of the Lorentz force and radiation pressure. We obtain its profile, 
both inside and outside the fast magnetosonic critical point, when the poloidal magnetic field is radial and 
\^ \ monopolar. Most of the energy flux is carried by the radiation field and the toroidal magnetic field that is 

wound up close to the rapidly rotating engine. Although the entrained matter carries little energy, it couples 
the radiation field to the magnetic field. Then the fast critical point is pushed inward from infinity and, above 
a critical radiation intensity, the outflow is accelerated mainly by radiation pressure. We identify a distinct 
observational signature of this hybrid outflow: a hardening of the radiation spectrum above the peak of the 
seed photon distribution, driven by bulk Compton scattering. The non-thermal spectrum - obtained by a Monte 

\ Carlo method - is most extended when the Lorentz force dominates the acceleration, and the seed photon beam 

is wider than the Lorentz cone of the MHD fluid. This effect is a generic feature of hot, magnetized outflows 
interacting with slower relativistic material. It may explain why some GRB spectra appear to peak at photon 
energies above the original Amati et al. scaling. A companion paper addresses the case of jet breakout, where 
diverging magnetic flux surfaces yield strong MHD acceleration over a wider range of Lorentz factor. 
Subject headings: MHD — plasmas — radiative transfer — scattering — gamma rays: stars 



. 1. INTRODUCTION 

Some astrophysical sources release energy at prodigious rates, as a result of their extreme magnetism and/or rotation. Gamma- 
ray bursts (GRBs) and soft gamma repeater flares are familiar examples, but the range of possibilities extends to rapidly rotating 
white dwarfs formed in binary mergers, or young and hot magnetars. The 'engine' generates an outflow that contains both an 
intense magnetic field and also a powerful flow of radiation . A gamma-ray burst jet beco mes hot as it works its way through a 
£Q ■ preceding layer of stellar material, as in the collapsar model (MacFadyen & Wooslev 1999); or interacts with a powerful neutron- 
^— 5 ' rich wind that was generated in the initial stages of a binary merger (e.g. Dessart et al. 2009 ). A blackbod y component is also 

\ expected to be carried outward from the engine itself even in the absence of jet interaction dGoodmanll986l) . 

We have previous ly argued that the full GRB pheno menon cannot be captured by 'thermal fireballs' dGoodmanl 119861 
emi & Piranl ll990). or 'Poynting-dominated outflows' (ILvutikov & Blandfordl 120031) alone. Instead both components are 
needed: thermal photons to provide seeds for g amma-ray emission, and large-scale magnetic fields to drive the outflow and then 
trigge r non-thermal activ ity and fast variability (lThompsonll994t Meszaros & Reesl[l997tlDrenkhahn & Spruit 2002; Thompson 
H ' 120061 [Giannios & Spruit! 120071 iZhang & Yanll201 lh . In some other circumstances, such as white dwarf merger remnants, the 
photon luminosity is bounded by the Eddington value, but the spindown luminosity can be vastly greater. (It is, of course, 
possible to find situations in which this inequality is reversed: the photon energy flux in giant magnetar flares is strong enough to 
pull some magnetic flux away from the star, but the total energy flux is probably dominated by photons.) 



1.1. Acceleration 

Our first interest here is in how such a 'hot electromagnetic outflow' is accelerated. We focus on stationary, axisymmetric flows 
in the ideal MHD limit, but add a radiation field with a prescribed source radius that can have arbitrary size with respect to the 
light cylinder. The inertia of the magnetic field dominates that of the matter to which it is tied, so that the outflow can achieve 
relativistic speeds. The magnetofluid is accelerated by two mechanisms: the Lorentz force (which operates in cold MHD winds); 
and scattering off the radiation field (which operates in thermal fireballs). 

In the optically thin regime, the radiation field is self-collimating, so that a relativistically moving frame can be defined in 
which the radiation exerts a vanishing net force on the matter. Radiation pressure dominates matter pressure, and with the 
possible exception of a small region close to the engine, the radiation temperature lies far below the rest energy of the advected 
particles. The outflow can therefor e be assumed cold outsid e the transparency surface, and scattering operates in the Thomson 
limit. In this first paper, we follow Goldreich & Julian (119701) in restricting the poloidal magnetic field to a monopolar geometry. 
This limits the efficiency of MHD acceleration (the fast critical point of the cold MHD fl ow sits at infinity). A mu ch stronger 
Lorentz force can arise through differential decollimation of magnetic flux surfaces (e.g. iTchekhovskov et all 2009). but even 
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then some parts of a cold MHD outflow may not have the requisite geometry. A companion paper considers jet geometries and 
allows for differential decollimation of the jet with respect to the radiation field (Russo & Tho mpson 2012 , paper II). 

The efficiency of radiative acceleration o utsid e the fast magnetoso nic point was noted by Thompson (2006) in the context 
of GRBs. Previous work by iLi et aTl {1992) and Bes kin et al.l ((2004) considered the interaction of a relati vistic MHD outflow 
with a quasi-isotropic radiation field. In a first approximation, the outflowing matter feels a net drag force. iBeskin et al.l (120041) 
focused on a slightly decollimating outflow with fast MHD acceleration, and analysed the influence of radiation drag in integral 
form through the changes imparted to the energy and angular momentum. They considered changes in the shape of the magnetic 
surfaces, and therefore in the MHD acceleration rate, imparted by the radiation field, but the first-order effect addressed in this 
paper - the outward acceleration due to the strong radial anisotropy of the radiation field - was absent in their calculations. 

We focus on time-independent outflows, because acceleration by radiation pressure becomes important before effects as sociated 
with the radial structure of the flow. The acceleration of a magnetized slab in planar geometry was studied numerically by 

iGranot et al.l d201 ll) . The slab is initially static; its mean Lorentz factor quickly reaches cr X J 3 , where oq is the initial magnetization, 
and then continues to grow. This effect may be relevant to GRBs after a jet has made a transition to a planar geometry at a distance 
r > cAf ~3x 10 n (Af /10 s) cm, where At is the duration of the prompt phase. If the shell is already moving relativistic ally, only 
the outermost fraction ~ 1 /2T 2 of the shell will experience this type of accel eration. The inte raction of this thin layer with an 
external plasma shell, which absorbs momentum from the magnetic field (e.g. Thompson 2006), must also be taken into account. 
Typically the bulk of the shell receives momentum from the radiation field before it comes into causal contact. 

We do not address how the outflow might accelerate below its scattering photosphere. For example, periodic reversals of a 
toroidal field might induce magnetic reconnection , break the degeneracy betw een magnetic pressure gradient and curvature forces, 
and thereby induce a net outward acceleration (Drenkhahn & Spruit 2002). The high efficiency claimed for this mechanism 
depends on neglecting the input of enthalpy, and therefore inertia, into the magnetofluid as the magnetic field reconnects. When 
an enthalpy source term is included, one only obtains relativistic bulk motion if a high fraction of the unsigned magnetic flux 
is removed by reconnection. This requires, at a minimum, that the net (signed) magnetic flux carried by the outflow is small 
compared with the unsigned flux, and that the non-radial magnetic field maintains a strictly uniform (e.g. toroidal) direction 
across current sheets. Our focus here and in paper II is therefore on ideal MHD effects. 

1.2. Non-thermal Spectrum 

The second focus of this paper is on the non-thermal spectrum of photons that are scattered by the outflowing matter. We 
ignore any effects of internal dissipation in the outflow, and focus purely on the spectral signature of bulk relativistic motion. We 
find a broad and flat extension of the seed spectrum when the outflow is rapidly accelerated by the Lorentz force, so that the seed 
photon beam is wider than the Lorentz cone of the MHD outflow. The smoothness of the scattered spectrum - in particular, the 
presence or absence of a residual bump at the seed thermal peak - is shown to depend on the net optical depth that is seen by the 
broader, unscattered photon beam. 

Existing calculations of non-thermal 'photospheric' emission from GRB outflows dGiannio s 2006; Beloborodov 2011) gener- 
ally assume that the Lorentz factor of the outflow has saturated at the t ransparency surface, as do calculations of the low-frequency 
spectral tail of thermal photospheres observed at oblique angles (e.g. lPe'eiil200a;lLazzati et alj|201 ll) . We broaden this approach 
by noting that locking between the advected thermal photons and the relativistic outflow will be broken in the presence of a 
somewhat slower component of the outflow. For example, material with a Lorentz factor ~ 3 - 5 times smaller than the mag- 
netofluid would be present as the result of the interaction of a jet with a star ( Thompson 2006). In a neutro n star binary merger, 
a relativistic magnetofluid could interact with a sub-relativistic, neutron-rich wind (Bucciantini et al. 2012). The slower material 
scatters advected X-ray photons into a broader beam, which continues to interact strongly with the faster magnetofluid even at 
low optical depth. It has long been realized that relativistic material movin g into an isotropic bath of very low-freque ncy (optical- 
UV) photons will upscatter them as it loses energy to Compton drag (e.g. Si kora et al.lll 994b iGhisellini et al.ll2000l) . In that case 
a spectral slope F v ~ i>~ 1 ! 2 is generated by the decaying peak energy of the upscattered photons - a different effect from that 
considered here, and a spectrum somewhat softer than that observed in the low-energy tails of GRBs. 

1.3. Plan 

The plan of the paper is as follows. Section [2] explains the efficiency of radiation-driven acceleration beyond the fast mag- 
netosonic point, and explains further our critique of acceleration by magnetic reconnection. The relativistic wind equations 
including radiation pressure are described in Section |3] and in Section|4]they are specialized to a monopole magnetic field. The 
flow properties near the fast critical point are analysed and the numerical method described. Numerical results are presented 
in Section [5] The spectrum of scattered photons is calculated in Section [6] by a Monte Carlo method, and the low- and high- 
frequency components of the spectrum discussed in the context of GRBs. Our results are summarized in Section|7] and the effect 
of rotation in the photon source is briefly discussed in Appendix lAl 

2. ACCELERATION OF MAGNETICALLY DOMINATED FLOWS BY RADIATION PRESSURE 

The outflow may be divided into an inner, optically thick part, in which the radiation field is effectively tied to the matter; and 
an outer part through which the radiation can flow almost unimpededQ Within this outer transparent zone, the radiation field 
becomes progressively more collimated with increasing distance from the engine. Even at the transparency radius, it can remain 
so intense as to push matter to much higher Lorentz factors than were achieved by hydromagnetic stresses operating in the inner 
zone. 

1 A h ard radiation spectrum, extending above m e c 2 in the frame of the ambient material, can trigger an e + e~ cascade outside the photosphere of a relativistic 
outflow I Thompson & Madau 2000; Beloborodov 2002). We ignore the effect of pair creation on the optical depth in this paper. 
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The radiation has a positive accelerating effect outside the scattering photosphere if the terminal Lorentz factor of the 
matter exceeds the Lorentz factor T(r T ) at the transparency radius r T . The intensity of the radiation field is measured by the 
compactness, which here is normalized to the mean mass m per scattering charge: 



xOv) = 



fnc^r T dfl 



(1) 



Here L 7 is the isotropic radiation luminosity, and <jj is the Thomson cross-section (we only consider classical electron scattering 
in this paper). Defining the scattering optical depth 



r,„ = / r^r^(l- (3 r )a T dr 
m 



(2) 



in a slowly rotating outflow with speed /3 r c and proper mass density p, the surface r es = 1 can be related to the flux of rest mass 
M, 

dM 



CJj 



6T 2 (r T )mc dfl 



(3) 



The magnetization of a relativistic outflow can be expressed in terms of the ratio of Poynting flux to rest energy flux, which is, 
far outside the speed of light cylinder, 



B: 



dM 



dLp "0 , in,, , 

~ — -r c ~ a c . 

dVL 4tt dft 



(4) 



The magnitude of a in GRB outflows is unknown, but there are a number of reasons to expect a range of values within dif- 
ferent components of the same outflow. The simplest case - but perhaps not the most relevant for classical GRBs - involves a 
iting i 

dMetzgeretal .1120111) . 

At r = r T , one therefore has the relation between photon compactness and magnetization, 



rapidly rotating nasc ent magnetar, whose wind can achieve a magnetization ~ 10 3 at late times after an early dirty wind phase 



The terminal Lorentz factor of a radiation-driven baryonic wind is (see Section |5^2| i 



roo(r T )~[x(rr)r(,v)] 1/4 = 



3 dL^/dCl 



6r J (r T ) 



dL P /dD, 

Hence the outflow emerges into transparency below this limiting Lorentz factor if 



1/4 



r(#v) 



r(r T ) 



6(L 7 /Lp)cr 



1/4 



<1, 



(5) 



(6) 



(7) 



corresponding to T(r T ) < a in a hot outflow with L 7 ~ Lp. Radiation plays a key role in accelerating the outflow when this bound 
is satisfied. 



2.1. Strong Radiative Acceleration Outside the Fast Critical Point 

An essential feature of a steady magnetohydrodynamic (MHD) flow is that the equation for the fluid speed, obtained by 
combining the continuity and momentum equations, becomes singular where it matches a normal mode of the fluid. In the 
absence of significant matter pressure, the relevant normal modes are the Alfven mode and the fast mode. 

In a relativistic outflow, the inertia of the advected toroidal magnetic field rapidly becomes insignificant outside the fast critical 
point. Here we consider the case where the poloidal magnetic flux surfaces do not experience the differential bending needed for 
efficient MHD acceleration. Then, when the radiative force is strong enough to provide significant acceleration at the fast point, it 
also controls the terminal acceleration of the outflow, as if the magnetic field were not present. In spite of this, the Poynting flux 
carried by the magnetic fiel d can continue to dominate the kinetic energy flux: the matter acts as a couple between the radiation 
field and the magnetic field (Thompson 200g). 

To see this, focus on the zone far outside the Alfven critical point. The solution to the induction equation is (Ferraro 1937) 

B<f, v^-n/rsing 

IT ~ " . W 

where 17/ is the rotational angular frequency of the magnetic footpoint at the boundary of the engine, and is constant along a 
magnetic flux surface. Equation dH) implies a uniform rate of transport of toroidal magnetic flux 



■5^. = v r rB^, ~ — CI sin 9r 2 B r = const, 



(9) 
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since the toroidal velocity has decreased to <C v r ~ c. We now allow an external (e.g. radiation) force to be applied to the fluid 
while imposing this constraint, as well as the constancy of mass flux 



dn 

The change in the radial energy flux is, at a fixed radius, 



dM , 

— = Ypv r r = const. (10) 



, 2 2 E e B^ \ 1 (dM 2 ST dM 2 % 



The right-hand-side here is singular at the fast point, where 

r 3 = P? = - ^ = cr. (12) 

c c 3 dM/dn 

The (toroidal) magnetic field dominates the inertia of the outflow in between the Alfven and fast points, but the small flux of 
ordinary matter dominates outside the fast point. 

Our main focus is on flows that are sub-fast magnetosonic at the transparency surface, and are accelerated through the critical 
point by the radiation force as well as the Lorentz force. For our adopted field geometry, this occurs when T(r T ) < er 1 / 3 , and 
requires a strongly magnetized outflow. 

A second type of flow is super-fast magnetosonic at the transparency surface, having already been accelerated through the 
critical point, T(r T ) > a 1 / 3 . The dominant acceleration mechanism outside the photosphere then depends on the poloidal field 
geometry. When the poloidal flux surfaces do not diverge from each other, the case studied in this paper, the Lorentz force 
freezes out and the remaining acceleration is by the radiation force. However, in a jet geometry there is a much more extended 
competition between MHD and radiation stresses outside the fast critical point, as we discuss in paper II. 

2.2. Insensitivity to Flow Structure Near the Alfven Critical Point 

Since the Alfven mode has the lower speed, its critical point sits closer to the rotating 'engine'. Inside the Alfven critical point, 
the magnetic field lines are effectively rigid and guide the outflow of matter and radiation. Outside, the rotation of the engine 
causes the magnetic field lines to be bent back into a Parker spiral, and the magnetic field is mainly toroidal. After a choice is 
made for the poloidal field profile, we wish to test the insensitivity of the large-distance flow solution to details close to the Alfven 
point. In fact, we obtain an essentially unique outer solution without explicitly requiring it to pass through the Alfven point. If 
the solution were not unique in this way, it would be suspect: the magnetic field configuration near the Alfven point maps directly 
onto the engine, and details of the inner flow solution depend, in turn, on details of plasma heating close to the engine. 

2.3. Bulk Acceleration by Magnetic Reconnection: A Critique 

The acceleration of a cold, spherical MHD outflow is inefficient due to a cancellation between the m agnetic pressure gradient 
and curvature forces. Rather than considering deviations from radial flow, Drenkhahn & Spruitl (120021) suggest that reconnection 
of a reversing toroidal magnetic field would convert Poynting flux to a large-scale relativistic expansion. Inside the scattering 
photosphere, they impose conservation of energy flux per steradian, 



dL 2^2 

— — = rrv, 

dn 



h+ (B,/T? 



4k 



(13) 



Here h ~ pc 2 is the bulk-frame enthalpy per unit volume, p is the proper rest-mass density, and B^/T the bulk-frame toroidal 
magnetic field. A sink term for the magnetic flux is imposed, cd r (rB^,v r ) = -rB^v r /r, but no corresponding source of h. With 
these assumptions, the dissipation of Poynting flux must be compensated by a growth in kinetic energy. A small decrease in 
Poynting flux corresponds to a large increase in T in a cold, high-cr outflow. 

More realistically, reconnection heats the particles, and in a high-a outflow there is a large increase in the particle inertia. Even 
allowing for rapid radiative cooling, the photons are trapped by the outflow, and 



5h ■ 



(B^/T) 2 



(14) 



Reconnection also creates a radial magnetic field, through the appearance of multiple magnetic X-points. The inertia of this small- 
scale magnetic field also dominates the particle rest mass. Taking both of these effects into account, one sees that dissipation of 
even half the magnetic energy can generate only mildly relativistic bulk motion. 

3. OUTFLOW EQUATIONS 

Given the technical challenge involved in adding the radiation force, we are forc ed to make simplifying as sumptions about 
the angular distribution of the poloidal magnetic flux. In this first paper we follow Goldre ich & Julian! (|T970) in allowing the 
magnetic field to bend in the toroidal but not the poloidal directions - i.e., the poloidal field is purely radial. This approximation 
is, in fact, increasingly well justified for very intense radiation fields, e.g., those which are strong enough to force open magnetic 
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field lines inside the speed of light cylinder. It should also hold with reasonable accuracy close to the magnetic equator. The 
solutions are then labeled by the fixed rest mass and radiation energy fluxes. 
To the steady Euler equation for a cold, relativistic MHD fluid, we add a term representing the radiation force, 



pTv ■ V(Tv) = — [(V • E)E + ( V x B) x B] + ^F rad . 

4n m 

Taking the dot product of equation (fl3T l with the poloidal magnetic field B p gives 



(15) 



Brdr+^dt 



IV --^ 
r 



B r d r (rTv^ + ^-dg(smerv^ 
smO 



4irTpr 



B 



B r d r (rB^+^-dg(sm8B^) 



sin 



1 

+ — 

m 



B r F™ d +BgFg ad ) 



The projection of the Coulomb force onto B^vanishes in a steady, axisymmetric MHD outflow with 



•B p =0. The, 



(16) 

= 0, since then E • B = 



-component of equation dT3l > is 



a e (Tv0 sin6»): 



1 



B r d r {rB^+-^-d e (sm9B^ 



4irTpr 

The toroidal and poloidal components of B are related through the flux-freezing condition ((5). 
The photon emission radius r s serves as a reference length, and the photon compactness is also measured at this radius: 



r 



n 



fx. 



Xs 



AixrJhc 



(17) 



(18) 



In a GRB outflow, the photosphere generally lies outside the light cylinder of the rotating engine, so we take lu = Vtjr s /c > 1 in 
our calculations. For example, an outflow launched by a millisecond engine which fills an opening angle 9 ~ 0.1 at a breakout 
radius of ~ 10 10 cm will have a physical width ttfrsinO/c ~ 100 times larger than the engine light cylinderQ This means that 
only a tiny fraction ~ 10~ 5 - 10~ 3 of the solid angle of a laminar jet will rotate rapidly enough that our approximation breaks 
down. (In practice, a modest amount of turbulence in the jet could mix this thin inner cone with the much wider, slowly rotating 
annulus surrounding it.) In practice, for computational ease, we consider intermediate values of w. In the case of an isolated 
millisecond magnetar one might indeed have w ~ 1. 
The magnetization is defined by 



_ B 2 n 2 f r 2 sin 2 8 

A-KTpVrC^ 

which is constant in the monopolar outflow considered hereQ The kinetic, Poynting and radiation luminosities, TMc 2 
are related by 



Ba 



a 
f' 



1 



Xs 



(19) 
Lp and L 7 , 

TMc 2 ' Txuj sin 9 B r T' TMc 2 6T 2 (x T )x T T 1 (20) 

where in the last equality we have have assumed that the flow is relativistic [see equation ©]. When radiation is absent, the 
energy and angular momentum per unit rest mass, 

t i a B^ . 
p, = 1 — - ; C = I x sin 

x<jjs\n9 B, 

are conserved along field lines. 

3.1. Relativistic Radiation Force 

The radiation emanates from a spherical static 'emission surface' of radius r s (x = 1). We adopt simplified spectral and angular 
distributions: the unscattered radiation is monochromatic, uniform at angles < 9 < ir/2 at the emission surface, and streams 
freely outward. At x > 1, the cone of the radiation field contracts and 



Ba,(J 



B r uj 2 x sin 9 



(21) 



l v = IqVq8(v-Vq) for 9 <9 S = sin 1 ^— ^ = i 



I v = for 9 > i 



(22) 



This accurately represents an isotropically emitting star that is surrounded by an optically thin wind. It still produces qualitatively 
correct results if the outflow is optica lly th ick near the engine, and experiences nearly linear growth of Lorentz factor with radius, 
driven by radiation pressure (Section [3.2) . 



2 Early claims of measurements of jet opening angles based on temporal breaks in GRB afterglow light curve (e.g. Rhoads 1999) have been revealed to be 
somewhat ambiguous based on more complicated behavior seen in the Swift data (e.g. Zhang et al. 2006). Nonetheless in a few cases detailed fits with direct 
hydrodynamical mo delling are consistent with opening angles ~ 0.1 (e.g. van Eerten et al. 2013), and a range of a few around this value seems likely (e.g. 
Panaitescu & Kumar 2002). 

3 This definition, following IMIcEel [\%% and IGoldreich & JuliarJ fi9"70l) . differs in terms CKy^/ClfrsraBf) to both a = (c 2 dM / 'dU)-' dL P / 'dtt, and a = 
(c^dM/dQT 1 ^, where $^ = v,-rsm0fB^ is the advection rate of toroidal flux. With this definition a is independent of radius if the poloidal field is restricted 
to be purely radial, but the last two definitions are non-constant at O(vw,/f2/rsin0/). 
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Because lu >> 1 typically, a non-rotating emission surface is a reasonable approximation. The effect of adding modest rotation 
to the photon source is addressed quantitatively in Appendix|A] and the corrections to the flow solution are shown to be small. 

The radiation interacts with electrons (and positrons) via Thomson scattering, and the radiative force is taken to be unperturbed 
by this interaction. The radiation force per scattering charge in a global inertial ('lab') frame is related to the force F lad in the 
flow rest frame by 



(23) 



where primes denote quantities in the frame co-moving with the fluid. Letting k be the unit wave vector of the incoming photon, 
we have the transformations 



ll=I„T 3 (l-/3-kf; dv' = dvT(\-j3-k); 



dSl'-- 



dCl 



r 2 (i-/3-£) 2 



jj_ k+{T-\)(p-W-pV 
V(l-(3-k) 



The radiation force in the fluid frame is given by 



F r 



ad/ - 



dp' oj 
~dt'~~ 



I'kfdtt'di/ = 



T z (l-f3-k 



fcr^+Ci-r -1 ) 



(3-k 



13-13 



dCl, 



and combining with equation ( 1231 gives 



l-P-k) k-/3T 2 (l-(3-k 



dtt. 



(24) 



(25) 



(26) 



The radiation force 
tions R, P, by 



one finds 



can be evaluated analytically by integrating over the emission surface. Defining dimensionless func- 

- 2 - 2 

F™ A = X s—R(r,r); Ff = X s—P{r,T), (27) 
l+2u 2 



R = - u r+^^+/3 r (2T 2 +v 2 ) Ji--+pJ -r 



1 



8 2 2 M r 

3 ~~~ 



'i+2r 2 +v 2 l 



1 1 



i-i 



i-4 

X 



3/2 



3/2' 



where the radial and non-radial four-velocities are 

(«,y) = (r/3,.,r/3 ). 

At large radius and Lorentz factor, expressions (EFi and (|29| i simplify to 

1 r2 . D ^ h 



R: 



4x 2 T 2 12X 6 ' 



i i r 2 

2x 2 \2T 2 + x 2 + 6x^ 



(28) 
(29) 

(30) 

(31) 



The first term in the radial force is due to photons propagating nearly parallel to the fluid; the second is the drag caused by photons 
which are aberrated into the anti-radial direction by the relativistic particle motion. The photon field seen by the particle is nearly 
isotropic in a frame were R = 0. At large \s and x not too large, the matter therefore tends toward the equilibrium Lorentz factor 



3'/ 4 - =3'/ 4 



(32) 



This frame only exist due to the extended nature of the source. 



3.2. Photon Distribution at a Displaced Photosphere 

In many cases, the outflow moves relativistically at its scattering photosphere, so that the radiation field is already collimated 
at the base of the transparent zone. Then the emission surface becomes a virtual one. Setting the photospheric Lorentz factor to 
the equilibrium value (l32l . the emission surface is pushed inward to a radius 

r,e ff =3'/ 4 -^-. (33) 
T(r T ) 

A radiation field that emerges from a relativistically moving photosphere is not cut off sharply at an angle 9 S = r S eff/r T = 
3 1 ' 4 /F(r T ), but has a somewhat smoother cutoff. Consider, instead, an intensity that is isotropic in the matter rest frame, I' v , = 
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Iq Vq 8(y' — Vq), and compare with the top-hat spectral distribution (l22t by normalizing to a fixed flux in the frame of the engine, 
F = 2ir J dfij dvl v (p). Then the frequency-integrated intensity in the lab frame is 

f tfvS I(9 = Q) 3 3 / 2 F ( 3>/ 4 \ 

as compared with / = F /tt^ 2 for 9 < 9 S . 

These two angular distributions yield the same photon force for a static charge at r 3> r?,eff an d, by construction, the same 
Lorentz factor for which the radiation force vanishes. The power with which an electron scatters photons in the bulk frame is 
nearly identical as well: P' = 4-ir J dv' u T l' v , = 3ajF /4T 2 for the isotropic bulk-frame intensity, and 

for the top-hat distribution. The spectrum of photons scattered by a cold electron flow is proportional to P' , and differs in 
normalization by ~ 20% in the two cases. 

More generally, if the outflow interacts with slower relativistic material - as discussed by Thompson (2006) in the context of 
GRBs - then the beam of photospheric photons is scattered into a cone of width > 1 /Y. This motivates our considering outflows 
with i) significant scattering optical depth at r = r s ; but ii) Y(r s ) > 1. Given the potential importance of such additional effects 
for GRBs, we stick with the simplest model a photon emission surface with uniform intensity, with the understanding that this 
a virtual surface in some circumstances. See Section |6]for calculations of the scattered photon spectrum and the application to 
GRBs. 




3.3. Acceleration of an Unmagnetized Outflow 

As a test of our formalism, we revisit the acceleration of cold matter by an intense radiation field, setting a = 0. Then equations 
( fT5l ), ( f2Tb and d28| i combine to give 



dY 
dx 



= Xs(r+ v -p); 



dx 



and simplify further to 



i r 2 



^£ D Xs_ | , 

dx ~ Xs ~ 4x 2 { Y 2 3x* ' ' 



dx 



1 1 

2x 2 T \ Y 2+ x 2 



Xsl 



(36) 



(37) 
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for large x, T. The matter Lorentz factor tracks r eq until the local compactness \ ~ Xs/F^x drops below ~ T (see Figure[T|i. The 
asymptotic Lorentz factor is 

r oo ~1.0S6 X l J 4 (38) 

as long as baryon loading is sufficiently low (77 = L^/Mc 2 ^> x'^ 4 ); otherwise r m ~ 77. 

The radiation force has a stronger effect when the flow remains optically thick out to a large distance from the engine, and 
develops at least a moderate Lorentz factor before the photons begin to stream freely. 

In this situation, the constraint L 7 < Lp would correspond to a limit on the photon compactness as measured at the photosphere, 



X( r r) 



— ^1 



4Trr T m e c 3 



-Xs: 



namely, 



X(rr)<6T 2 (r T )a 



(39) 



(40) 



[see equation ([2011, Now the parameter \s describes a 'virtual' photon source that is buried at large scattering depth, and is 
already collimated at the photosphere. The terminal Lorentz factor increases to 



xl /4 = 



3" 1/4 r(r T ) X (r T ) 



1/4 



(41) 



If the fast point were to sit in the transparent zone, then T(r T ) < a 1 1 3 and one has the bound 

Xr <6a 5 / 3 ; r'oo.x ~ L3<7 1//2 . (42) 

If the bound ( f42l > is not satisfied, then the outflow will pass through the fast critical point before reaching its photosphere. 
Nonetheless, the radiation field will provide significant supplemental acceleration outside the photosphere for any value of T(r T ) 
below a. 

4. DIMENSIONLESS OUTFLOW EQUATIONS: 6 < vr/2 

The poloidal magnetic field lines make a transition from dipolar to nearly radial at the light c ylinder of an isolated pulsar, 
being forced open by causality constraints (IConto poulos 2005; McKinnev 2006; Soitkovskv 2006). More collimated outflows 
are expected from magnetized stars surrounded by accretion disks {Blandford & Znaiek 1977; Blandford & Pavn el 19821) . and are 
the focus of paper II. 

We consider a streamline close to the magnetic equator (but not so close as to involve the equatorial current sheet). We work 
with a simplified form of the Euler equations that includes the radiation force, but assumes the poloidal magnetic field to be 
purely radial. Therefore we substitute 



Br 



Be v 



B r ~^; -^ = -«l; 



sin 9 ~ 1 , 



into equations ( fT6l l, (fTTI i, which gives 



d r Tc z 



v <t> 



drirTvt) = --^hdArB^+^F^; 
Attl pr m 



B r 



-d r (rT V<p ) = —^dArB^-Fl 
r 47rl pr m y 



(43) 



(44) 



(45) 



The various terms in the right-hand sides of these equations can be separated into purely magnetocentrifugal pieces (which do 
not depend on the radiation force), the direct radiation force, and a cross term: 



dT 

dx 



r" +r +r' 

Meff 



dx 



where 



r' +r' +r' ~ 



f3 r xuj 



(3 r B r 



Meff 
°Xs 



1 + - 



bBj 



B r B r 



R+^P 

B r 



(46) 



(47) 
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X 2 0J 2 T 2 U 2 B r \ B r 



(48) 



and 



a , of v 2 \ 2ovT 
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4.1. Flow Through the Fast Critical Point: 9 < ir/2, B r ~ r~ 2 

These MHD wind equations have critical points where the flow matches the speed of a cold MHD mode. The Alfven critical 
point corresponds to a flow speed u = er(l - lo 2 x 2 ) / ui 2 x 2 , and therefore sits inside the light cylinder. 

Since we are looking for robust flow solutions which do not depend on the magnetic field structure in this inner zone (Section 

I2.2K our focus is on the fast critical point, where 

rr f - - R 2 \ 

(50) 

This is the singularity appearing in equation (I46t . A regular flow solution passing through this point must satisfy two conditions: 

T' a (xc) + T' x (x c ) + T' ax (x c ) = 0; Meff(* c )=0. (51) 

Here x c is the (so-far undetermined) radius of the fast point. Equations dBTT l generate a one parameter family of flow solutions. 

The specific angular momentum evolves according to the simple equation dC/dx = (x/(3 r )x s P- We have tested equations 
d44l) and d45l l by combining them and re-deriving this equation. This means that the critical point of the equation for dfi^/dx 
does not impose independent constraints on the flow solution: a solution with a smooth T(x) profile automatically satisfies the 
^,M) + ^ x (x c ) + ^ r7X (x c ) = 0. 

The presence of a cross term in the wind equations, involving both the magnetization a and the compactness Xs, deserves some 
comment. Inside the fast magnetosonic point, and outside the Alfven point, the radial motion of the fluid has effectively a negative 
inertia: a positive external radial force extracts energy from the outflow. This negative inertia arises from the response of B<j» 
which dominates the energy integral (|2TT >. to changes in the radial flow speed. The scaling B^ ~ f3~ l implies a decrease in toroidal 
field energy with increasing f3 r . On the other hand, the angular momentum (|2"TT i is also dominated by the electromagnetic field. A 
change in B^ creates unbalanced toroidal stresses, which are a source for fy, and are proportional to a\ s - Changes in f3<f, and /3 r 
of opposing signs allow B^ to remain nearly constant. The net change in Lorentz factor is positive, 5T = T~ 3 (f3 r S(3 r + (3^8 (3^) > 0, 
if 1 + P^B^/ f3 r B r > 0. Close to the light cylinder, where the term ~ -o/ux 2 lo 2 in p^ dominates the term -cr/w 3 , the effective 
inertia has the usual sign. 

4.2. Pure MHD Wind (Xs = 0) with Monopolar Radial Magnetic Field 

Whe n radiation fields are absent and r> 1 our equations reduce to the cold limit of the system studied by Goldrei ch & Julian! 
(I1970I) . The angular momentum C [equation (ETT ll is conserved, and determines 

xu>{Cu>l3 r - a) 



x uj u — a 

The fast point lies at infinite radius, the Lorentz factor being limited to its critical value 



(52) 



r 00 , (T = Vl + (T 2 / 3 ~a 1/3 . (53) 

The unique solution passing through the Alfven and fast critical points is the minimum-energy solution found by lMichel |l969). 
The slow acceleration of radial flows is an artifact of the near perfect cancellation of the outward magnetic pressure gradient 
force and the inward c urvature force. Faster acceleration is possible through a faster-than-spherical d ivergence of the outflow 
(Begelman & Li 1994), or differential bending of the poloidal field lines (Tchekhovs kov et all 12009). the effect of which we 
examine in paper II. 

4.3. Small Compactness Limit: Xs 4<t 4//3 u;~ 2 

For small but finite Xs, the Lorentz factor at the fast point remains unchanged from the pure MHD solution, T c ~ u 1 / 3 , but the 
critical point is brought in to a finite radius. The pure magnetocentrifugal and radiation terms in the wind equations dominate at 
large x, T, and simplify to 

K^~(2-^)^-^; r^ Xs R^j%- 2 (54) 

UX Z U> \ Pr J UX Z UJ * 4x z I L 

Applying the regularity condition Y' a + T' ~ at the critical point gives 

4cr 5 / 3 a 1 XsU - - ~ ° cr 1/3 



£ c ~r c x c /^. c +-J-~ — (l + a 2 ") (1ow Xj ). (55) 



C Xs^ 2 ' ' x c UJ 4(T 5 / 3 ' c c qi,c w 

4.4. Large Compactness Limit: Xs 3> 4cr 4 / 3 w -2 
For very large Xs the fluid is locked to the radiation field while crossing the fast point at a finite radius, 

r f ~r eq => x c ~3-'/ 4 r c (hi g h Xj ). (56) 

Here B r need not be small compared with B^. Approximating <C xuj and making use of equation dH), equation ( f50b becomes 
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Substituting u c ~ T c ~ 3 1,/4 x f gives the critical point parameters 



rV3 



1 + 



1/3 



i + -r- 2 

2 c 



(high x*)- 



(58) 



The transition between the low- and high-Xs scalings for the fast point occurs where jc c i ow ~ x c h; g h, at a compactness % s ~ 

4.5. Numerical Methods and Boundary Values 

The singularity at the fast point is difficult to handle with standard integration methods, especially when the wind equations 
become stiff, as they do when % s ^ 1. We employ two different integration schemes: a shooting method for low Xs> an d a 
relaxation method for high \s- The position x c of the critical point is unknown a priori, since the regularity condition provides 
only one constraint on the two coupled ODEs. Once a candidate value of x c is chosen, the flow variables at the critical point are 
uniquely determined, as is the flow solution interior to to it. However, this interior solution generally diverges at small radius. 
A strong divergence is avoided only for a narrow range of x c , and even then the solution tends to develop sharp gradients in T 
and /3$. Avoiding such gradients leads to an essentially unique choice of x c and a robust flow solution. We have found that these 
smooth solutions have the property that the terms in dfi^jdx which depend explicitly on the radiation field should sum to zero at 
x = 1 

This recovers Michel's minimum-energy solution (Michel 1969) in the limit of vanishing radiation field, Xs 0- The resulting 
initial values of T and are plotted in Figure |2]for a = 10 3 and uj = 2. 

w=10, cj=1000, j8 0R =O 

■ 1 1 ii^ 



T(x=i) 




0.3 r 
0.2 E- 
o.i F- 



mi^ 1 1 1 1 1 i limn 1 1 mm i ■ ■ i ■■■■[ 1 1 ■ i ■ ni| i Minn 1 1 1 1 1 m rrniBj 



^. eq (x=l) 3 



nmJ ■! J I I I iii I i i limit 



0.1 1 10 1 00 1 3 1 4 1 s 10 6 10 7 



FIG. 2. — Lorentz factor and rotation speed at the inner boundary of an o utflo w with a = 10 3 , ui = Qfr s /c = 2. Dotted line corresponds to a vanishing net 
toroidal acceleration from the radiation field at the inner boundary, equation 1591 . 



At low Xsj the flow solution interior to x c is obtained by first determining the locus of critical points, and the corresponding 
values of T(x c ), 0d> (x c ). We then sh oot inward from x = x c using a 5th-order Runge-Kutta algorithm with adaptive step size (see 
Sections 7.3, 7.5 ofQCiusalaas 2010). The value of C c - which, unlike x c , is single valued - is iterated until the required small-* 
behavior is obtained. This method fails when Xs 100, since the equations become extremely stiff near x = 1, and machine 
precision becomes inadequate to distinguish values of | x c ,r(x c ),/3rf,(x c )} tha t lead t o converging and diverging solutions. 

At high Xs, we use the relaxation method described in lLondon & Flannervl JT982) for transonic hydrodynamic flows. The wind 
equations are replaced by finite-difference equations on a grid. Starting with a simple trial solution, and an initial guess for x c , 
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the inner boundary condition ((59) is applied along with the regularity condition ( BTT i at the critical point. Since the location of the 
critical point is unknown, an independent variable q which labels mesh points is introduced, along with a mesh-spacing function 
Q(x) = (where is an unknown number). We choose Q(x) oc lnx, which tightly packs the grid points near x = 1 (where 
the equations are stiff), and spreads them out near the critical point (so as to avoid divergences induced by the singularity). 
This necessitates adding two more ODEs, for x(q) and ^. The error in the initial guess at grid point i is quantified in terms 
of Ej = v; - y,-i - (xi - Xj-t )dy/dx and a new solution is obtained by a multi-dimensional Newton-Raphson m ethod (see sectio n 
18.3 of Pre ss et al71 f2007). The matrix inversion is performed using a Gaussian pivoting algorithm detailed in iKiusalaasi (2010). 
The critical point must be approached from below with this method and so the trial solution is typically cut off at ~ 90% of the 
expected critical radius. This method fails in the low Xs regime since the critical point is at large radius where dT/dx ~ 0, and 
even a slight overshoot in T c will be catastrophic. 
In both the low and high x.i regimes, the solution is completed by shooting outward from the critical point. 



gj=10, (7=1000, /ff H =o 




10 
10' 



i i i i i 



10 

x = r/r 



100 



FIG . 3 . — Acceleration of a magnetized outflow (cr = 10 3 , oj = fi fr s /c= 10) for varying compactness of a radiation source situated at x = 1 : \ s = 10 2 — 10 7 
bottom-top on the right side in T(x), top-bottom in (S^(x). 



5. RESULTS 

As the initial radiation compactness Xs is increased, the plasma accelerates faster at small radius, and reaches a higher terminal 
Lorentz factor. The radiation field is weakly collimated near the emission radius, and so the starting Lorentz factor is reduced at 
large Xs- Figure [3] shows the dependence of V and on radius, obtained for a strongly magnetized outflow (cr = 10 3 ), starting 
from 10 times the light cylinder radius (u = 10). The Lorentz factor profile shows a smooth transition from Michel's minimum- 
energy solution in the low-x. s limit, to the purely radiation-driven solution at large Xs- An explicit comparison of magnetized and 
unmagnetized flows is made in Figure |4] Once sees magnetization at a level a = 10 3 makes little difference to T{x) for Xs ^ 10 5 . 

The reducti on in the starting Lorentz factor by radiation drag is tied to the negative inertia of the plasma inside the fast critical 
point (Section HTTT l. In this inner zone, the inertia is dominated by the magnetic field. A negative radial force, as is provided by 
the radiation field when V > r eq , pushes the flow to higher speeds. The situation reverses at the fast point, where matter begins 
to dominate the inertia and V drops below T eq . Outside the fast point, radiation provides a positive push on the matter and the 
entrained mag netic field. The asymptotic Lorentz factor grows with respect to the cold MHD flow, as we discuss in more detail 
in Section |5~3| 

The fast critical point (marked by the red cross in Figure |3} sits at infinity in the cold MHD wind with monopolar magnetic 
field, but moves in rapidly as Xs increases above 10-100. The explicit dependence of the flow properties at the critical point on Xs 
is shown in Figure|5]for a = 10 3 , and two launching radii (ui = 10,2). Expression d55l l provides an excellent fit to the fast critical 
radius x c for Xs ^ 10 3 , above which it settles to the value where the equilibrium Lorentz factor r eq = cr 1 / 3 ; see equation d32l . The 
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radial field inertia causes a small upward adjustment in T c , as may be seen by comparing the cold MHD expression T c ~ cr 1 / 3 
with equation (l58l l. There is a small downward adjustment in the total angular momentum at the critical point. 

The coupling of the radiation field to the magnetic field has some subtle effects on the rotation of the outflow. The angular 
velocity at the critical point decreases with increasing \s at large compactness, as would be expected from the increasing friction 
imparted by radiation field. (The radiation field is assumed not to rotate; the effect of rotation is considered in Appendix lAl) 
However, /3$ iC increases with \s at small compactness, due to the shrinkage in the critical point radius. The finite value of /3<j, iC is 
also worth commenting on. The outward acceleration of the flow requires a finite angular speed for the matter. Strong radiation 
drag, acting alone, would rapidly damp the matter rotation. However, the radiation field enters indirectly into equation (l48l l for 
dp,/, jdx < through the term proportional to a\ s . This cross term represents the reaction of the radiation on the Lorentz force. It 
is generally positive and almost exactly cancels the azimuthal radiation drag. (We enforce this condition at x = 1, but the condition 
is nearly satisfied automatically at all radii.) This near cancellation allows the matter to maintain a high enough azimuthal speed 
to reach the fast critical point. A further interesting effect, evident in Figure [3] involves the near constancy of (3^ at small radius 
in radiation-dominated outflows (\ s 3> cr). This is due to a near cancellation of all the terms in dfitp/dx. 

5.1. Exchange of Energy Between Radiation and Magnetofluid 

In a hot magnetized outflow, the fluid acts to couple the radiation and magnetic field allowing energy to be exchanged between 
them as the flow accelerates. To study this exchange we write the total energy in terms of kinetic, Poynting and radiation 
luminosities as 

Lk + Lp + L^ 1 B xx(x) 

Mc 2 xuj sin o B r OT eS)& l L s 

where we have taken V oc r in calculating the optical depth. For simplicity, our solutions (being calculated in the optically thin 
regime) assume a constant radiation flux and thus do not strictly conserve energy. Nevertheless, we can impose conservation 
using the obtained profiles for the kinetic and Poynting fluxes to study the qualitative features of the energy exchange. The results 
are shown in Figure[6]for various values of photon compactness while taking the optical depth at r = r s to be unity. The radiation 
luminosity generally increases sharply at small radius where photons are upscattered by the highly relativistic flow. This must 
be accompanied by a corresponding decrease in the magnetic energy since the Lorentz force is still accelerating the flow. To see 
why this occurs note that well outside the light cylinder the Poynting luminosity in a monopolar outflow is Lp ~ f3~ l aMc 1 which 
can decrease significantly only if the flow is not yet extremely relativistic at the photosphere. In this case field lines are coiled 
tight and transfer their magnetic energy to the radiation field as they unwind. The matter acts mainly as a catalyst in this process, 
with relatively small changes in kinetic energy. 

Outside the fast point the Poynting flux is essentially constant and further changes in the kinetic energy come directly from the 
interaction with the photon field. (In Figure |6]the positive/negative slopes of the energy profiles are plotted as dotted/dashed the 
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FIG. 5. — Top panel: Radius x c of the fast critical point, and corresponding flow variables V, fi^ and C, as a function of radiation compactness \s, for tr = 10 3 , 
ui = 10. Red and blue dotted lines show the analytical approximations 455 1 (low \s) and 45 8t (high \s)- Bottom panel: Now with a smaller launching radius, 
ui = 2. 
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FIG. 6. — Energy exchange in an outflow with a = 1000, uj = 10 with an optical depth of unity at r = r s . Poynting, kinetic and radiation luminosities are plotted 
as solid lines and dL/dlnr for each is shown in corresponding colors with dotted/dashed representing positive/negative slopes. 



corresponding colors.) When \ is high* the fast point marks a sharp transition in the energy exchange: conversion of magnetic 
to radiation energy inside; conversion of radiation to kinetic energy outside. Accounting for the increase in photon luminosity at 
small radius would enhance the acceleration outside the fast point. 

5.2. Application to Gamma- ray Burst Outflows with Displaced Photospheres 

We briefly consider the application of our results to gamma-ray bursts. When a jet emerges from a Wol f-Rayet star, or from 
a clou d of neutron-rich debris, much of the jet material may experience a strong outward Lorentz force (iTchekhovskov et all 
2009). Even in that circumstance, the radiation pressure force can enhance or impede the acceleration of the jet material outside 
its photosphere (Russo & Thompson 2012, Paper II). Other parts of the jet may have unfavorably curved flux surfaces and feel a 
weaker Lorentz force, more typical of the poloidal field geometry examined in this paper. 

Although our solutions are found in the zone exterior to a static, photon-emitting surface, the photon field that emerges from 
a displa ced photosphere will have a similar effect on the outflow exterior to it, if we define an effective source radius r sfi s as in 
Section 13.21 A high photon intensity also has the effect of driving the flow profile to a linear relation T(r) oc r inside the fast 
critical point. In that case, the flow profile does not depend on whether the outfl ow is optically thick or thin inside the critical 
point. The conditions for this to be the case are outlined at the end of Section [3~3l 

To fix some numbers, consider a gamma-ray outflow with an (isotropic) luminosity AirdL^/dVL = 10 51 L5i erg s" 1 and a photo- 
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sphere of radius r T = 10 I0 r r ,io cm. The corresponding photospheric compactness <[T} is x( r r) = 1 x lQ & L5ir T l m (m p /m). Trans- 
parency at r < r r is guaranteed if T(r T ) > 10, and pairs have mostly annihilated. The effective source radius, obtained by equating 
the bulk Lorentz factor (l3Zt of the radiation field with T(r T ), is r s , e ff ~ r T /T(r T ) < 10 9 r r ,io cm. The compactness scaled to this ra- 
dius, which determines the amplitude of the radiation force in the wind equations d46ll-(|49l), is Xs ^ L5ir~ l m [T(r T ) / 10]{m p /m). 

5.3. Asymptotic Lorentz Factor: Fully Transparent Outflows 

Since we are considering a strictly monopolar poloidal magnetic field, the asymptotic Lorentz factor is limited to ~ a 1 / 3 
at small radiation compactness. Higher Lorentz factors are possible at large Xs, where acceleration is dominated by radiation 
pressure. First consider the case where the outflow is fully transparent at r ~ r s , and moves transrelativistically near the inner 

boundary. Then T oox ~ 1.086x' ,/4 (Section 13.31 ). In the example at the end of the preceding Section, this corresponds to 
1\ A - 200. 

We find a smooth transition between these two limits as shown in Figure|7]for a = 10 3 and w = 2, 10. This is well described by 
the function 

^oo — -^oo,ct ^^oo.x' (61) 

with n ~ 2.5. In this situation, strong radiative acceleration requires a photon source that is not sourced internally by the outflow. 
Otherwise, the photon luminosity is approximately bounded above by the Poynting luminosity, which implies L 7 < Lp and Xs < cr. 




0.1 1 10 10 2 10 3 10 4 10 5 10 6 10 7 1 10 10 2 10 3 10 4 10 5 10 6 

FIG. 7. — Asymptotic Lorentz factor of a magnetized outflow (<r = 10 3 ) exposed to a central radiation field with compactness \ s at x = 1- Left panel: ui = 10. 
Right panel: 10 = 2 



6. SPECTRUM OF SCATTERED PHOTONS 

The outflowing matter scatters the radiation field, and the frequency distribution of scattered photons is non-thermal. The effect 
grows stronger as the seed photon beam grows wider with respect to the Lorentz cone, 6 S > 1 /T, We now calculate the spectrum 
in a situation where most photons see a low optical depth after scattering by the magnetofluid. 

To motivate our calculation, it is worth discussing how such a situation could arise in the context of GRBs. P hotons which 
start with a nearly blackbody spectrum will maintain that spectral distribution in a simple relativistic fireball (e.g. Belo borodovl 
1201 UlLazzati et al]|201 lb . At the photosphere, the photon beam has an opening angle 9 S ~ T -1 and maintains a nearly isotropic 
distribution in the bulk frame. So we expect minimal frequency redistribution in MHD outflows that are accelerated mainly by 
radiation pressure, V ~ r eq . 

The outflow feels a very strong Lorentz force where neighboring magnetic flux surfaces diverge from each other, and we show 
in Paper II that this can driven V > r eq even close to the photosphere of a relativistic outflow. But even if this effect is not 
operating, there is still good reason to expect that material from the progenitor star (in the case of collapsars), or from a preceding 
neutron-rich wind (in the case of binary NS mergers) will interact with a relativistic MHD outflow. This material can be massive 
enough to maintain a Lorentz factor r s i ow much lower than that of the MHD outflow, and so broaden the photon beam into a cone 
of width 9 S ~ 1 /r s iow F° r example, a precursor shell, trapped at the head of a relativistic jet, becomes Rayleigh-Taylor unstable 
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as it is pushed outward, opening out zones of an angular size < 1/T in which the jet material can flow freely (Thompson 2006). 
Another possibility is a jet boundary layer containing material of an intermediate Lorentz factor. 

It is, therefore, posssible to probe the spectral signature of slow material in a GRB, independent of the details of photon 
creation at large optical depth. Indeed, because photons moving off the axis of a relativistic flow scatter at a higher rate, this type 
of interaction can occur even outside the photosphere of the MHD outflow, and still result in significant rescattering by the faster 
material. There are interesting implications for both the high- and low-frequency portions of the photon spectrum. 

6.1. Scattering of a Monochromatic Source with Uniform Intensity 

We are interested in the power radiated from a steady flow, with a fixed particle density at a given radius. In the case of an 
impulsive event such as a GRB, of duration At, the flow is effectively steady near the transparency surface if T 2 (r T ) ^> r T /cAt. In 
this situation, there is no correction between the time coordinate of an observer sitting at a large (but non-cosmological) distance 
from the engine, and the time coordinate t of the engine rest-frame. 

We start with a monochromatic photon source I u = Iqvq8{v-vq), and the top-hat angular distribution (l2Zt corresponding to 
a (virtual) emission radius r s . More general source spectra are then considered by a convolution. The condition of elastic 
(Thomson) scattering is 

The direction cosines of the incident and emitted photons are measured with respect to the radial direction (we neglect any small 
non-radial motion of the matter), fi = k-r, /x em = k em ■ r. They take the range 



1 



<fi<l; 



"I < Mem < 1- 



(63) 



Primed quantities are measured in the rest frame of a scattering charge. 
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FIG. 8. — Maximum frequency of scattered photons (upper curves) and minimum frequency (lower curves), as a function of radius, in a strongly magnetized 
spherical wind (a = 10 3 ). Compactness of photon source at x= 1: Xs = 1 (solid black); Xs = 10 2 (dotted blue); Xs = 10 4 (short-dashed green); Xs = 10 6 
(long-dashed red). 



A scattered photon reaches the maximum frequency 



l-/3Mmin A , T 



^2 



.,,„,_ ( _ ? -v ~ ( v 1 + ^2"J ^ (^ r >!)- (64) 

Hence the scattered spectrum develops a significant tail at frequencies above uq if the outflow has a Lorentz factor T > r eq 
[equation (l32ll. This non-thermal tail is more pronounced at lower values of the compactness x.s, where the Lorentz force 
dominates the acceleration of the flow. The scattered spectrum also extends to a low frequency tWinm = [(1 ~ /?)/(! + PfWo ~ 
vq/AT 2 . Sample profiles of t'em.max and v em ^ m are given in Figure|8] and the asymmetry of the rest frame photon flux in Figure|9] 
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FIG. 9. — Asymmetry of the target photon distribution in the bulk frame of the outflow, represented by the minimum direction cosine of th e photons //, as a 
function of radius, pi is measured with respect to the flow direction, and is larger than the plotted value at very low frequencies [equation A77H . Compactness of 
photon source at x = 1: \s = 1 (solid black); Xs = 10 2 (dotted blue); Xs = 10 4 (short-dashed green); \ s = 10 6 (long-dashed red). 

The frequency v im of the emitted photon is hardest in the flow direction. Softer photons are emitted off this axis, and also see 
a larger optical depth to scattering. In a spherically symmetric outflow, a photon emitted at radius r sees a scattering depth 



o T dN [°° dr 2 

7" es (r,/W = — — / [l-p(r 2 )Wz)] at w \ 2 - 
c dil J r p(r2)fJ.(r2)r2 Z 

The direction cosine evolves from the emission radius r to r2 > r according to 

2 



(l-Mem) (r 2 >r). 



(65) 



(66) 



It is straightforward to calculate the emergent photon spectrum by a Monte Carlo method. The direction cosines of the input 
photons are drawn randomly from the uniform distribution d63l l at x = 1. The optical depth Ar es to the first (next) scattering is 
determined by randomly picking 1 -e" ATcs , followed by a step-by-step integration of r es along the ray. Scattering is performed 
in the rest frame of the cold flow, by transforming // = (fi — j3)/(l — picking rest frame scattering angles 8' s , (j>' s with respect 
to this axis, and then determining the direction cosine of the outgoing photon via fi' em = fj/ cos6>' s + (l- /i'^'^sinf^cos^. The 
photon escapes if Ar es exceeds the total optical depth along the ray. Working in spherical symmetry, we record the frequency but 
not the direction of the outgoing photon. 

6.2. Linear Flow Profile 

To illustrate some of the main effects, we show in Figure [K)]the spectrum resulting from an outflow with a simple linear profile 
T(r) = 5(r/r s ), and various values of the scattering depth experienced by the most obliquely propagating photons. In this case, 
the bulk frame of the seed photons moves at a somewhat lower Lorentz factor than the magnetofluid, r s i ow ~ 0.2 V. The scattered 
spectrum therefore peaks well above the seed frequency, ^ e m,max ~ 30 uq [equation d64l)l, driven by bulk Comptonziation. 

First consider the output from a monochromatic source spectrum (the black lines in Figure [Tot. The output spectrum is fairly 



flat, F v 



,0.5 



, over a decade in frequency below the peak, steepening to F v ~ v at v ~ i/q. The scattered photons see a small 



optical depth (suppressed by a factor ~ (r s i 0W /T) down to fairly low frequencies, where the slope approaches a Rayleigh-Jeans 
value. 

As the optical depth of the seed photons is increased, the spectrum steepens slightly at high frequencies. The spectra in Figure 
( TTOb are labelled by the optical depth seen by seed photons at the maximum angle 9 S ~ (r/r s )~ l and minimum direction cosine 
/i m ; n ~ 1 — 6^/2. Requiring that T es (r s , /z m ; n ) > 1 allows the radial optical depth to remain small if T ■ 6 S ^> 1: 

T e s(r s ,fJ,min) ~ [l + r 2 (r s )l T es (r s ,l) T es (r s , 1) > A ■ (67) 
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env env 

FIG. 10. — S pect rum of scattered photons emerging from a spherical outflow with a simple profile Y(r) = 5(r/r s ), obtained by direct Monte Carlo integration. 
Optical depth l !67t of the off-axis photons is 0.3, 1, 3, 1 0, 30 from bottom to top, corresponding to optical depths < 1 after scattering at frequencies > uq. 
Unscattered source photons not included (see Figure HTI for comparison). Black lines: monochromatic source spectrum. Left panel: black-body photon source 
(green lines). Right panel: GRB-like source spectrum (red lines), F u ~ const X e~ h "l kTa (the low-frequency half of the Band function extended to higher 
frequencies). In both panels, the source spectrum is the dotted curve. 
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FIG. 12. — Total spectrum of photons emerging from a hot electromagnetic outflow (tr = 10 3 ) with a blackbody source of compactness Xs = l" 2 (left panel) 
and Xs - 10 6 (right panel). Flow profile shown in Figure [3] Curves from bottom to top have off-axis optical depth fl67t equal to 0.3, 1,3, 10,30. Increasing 
i"es('"j,A t min) implies that the outflow moves closer to the equilibrium value r cq at the transparency radius, hence the steeper high-frequency spectrum. Dashed 
spectra do not include unscattered source radiation. 

6.3. Generalization to Blackbody and Other Thermal Seeds 
It is straightforward to calculate the output spectrum for a broadband source with spectrum F v q{vq) via 

' v e . 



em 



dv' Q 



(68) 



Here F Vem (i> em ) is the Greens function response (|79l of the scattering outflow to a line photon source. 

The result for a blackbody seed is shown in the left panel of Figure[l0] and for a thermal seed with a GRB-like spectrum at low 
frequencies, F v = constant x e - lm / kT <>^ j n the right panelQ In both cases, the seed temperature is normalized so that F v peaks at 



1 This corresponds to the low-frequency half of the Band function I Band et al. 1993). The high-frequency power-law tail of the observed GRB spectrum must 



20 



Russo & Thompson 



c 



10 2 
10 

1 



J 0.1 



0.01 

io- 3 



10 



-4 



E 1 1 Mill 


i i nun i i nun i i mini 


1 Mill 




X s = 1-10 6 (solid) 




=- 


T(x) = 5x (dashed) 






T es( r s'M min ) = 10 




— 


v^^^^xX \ \ * 

'■ \\ \ \ 1 
•tir \\ \ \ v 
'ir \\ \ \ ' 


llllll 1 1 Mill 1 1 1 





1 ; Jil in 11 


1 1 llllll 1 1 = 



0.01 0.1 1 10 10 2 10 3 



em' 



V 







FIG. 13. — Same as Figure [T2l with a GRB-like photon source (dotted curve), for a range of compactness \ s = 1, 10 2 , 10 4 , 10 6 (solid curves, softening at high 
frequencies with increasing compactness). Dashed line shows result for linear flow profile, T(x) = 5x. Outflow has a fixed off-axis optical depth T es (r,,/i m j n ) = 10. 
See Figure[3]for corresponding flow profiles. 

v = vq. In the second case, the upscattered portion of the spectrum at v > vq connects smoothly with the seed spectrum at v < vq. 
One can also add in the source spectrum, appropriately attenuated, 



F vQ -> 2tt / I„ Q exp[-T es (r s , n)]dn 



i+r- 2 (r s ) 



I exp [-r(r s , 1 )] - exp [-r(r s , fi min )] j F u0 , 



(69) 



7"es(ry; A^min) 

where we have made use of the small-angle approximation (l72t to the scattering depth. The effect is shown in FigureQT] The seed 
blackbody peak is apparent if the seed photons see a maximum optical depth r(r s , /j, m i a ) ~ 1, and helps to extend the low-energy 
tail for r(r s , //min) ~ 3. It is subdominant for larger optical depths. 

6.4. Results for Radiatively Driven MHD Winds 

In Figures [T2lfT5l we show spectra calculated using the outflow profiles obtained in Section|5] with flow parameters a = 1000, 
uj =10. A distinct feature of the spectra is a prominent high energy tail appearing in outflows with moderate \s ^ 10 3 - compare 
the spectra for \s = 10 2 and 10 6 in Figure [12] At low compactness, the outflow acceleration is dominated by MHD forces, and 
T 3> r eq at the base of the outflow. When the acceleration is dominated by radiation pressure, increasing the optical depth only 
causes small changes in the peak of the scattered spectrum (as for Xs = 10 6 ). A direct comparison of outflows with a range of Xs 
is made in Figure [13] 

The spectral index of the Xs - 10 2 outflow is shown in Figure [14] for both blackbody and GRB-like seed spectra. As the 
optical depth at the base of the outflow increases, the cutoff frequency drops and the spectrum softens. In the blackbody case, 
the low -ene rgy tail has a fairly constant, Rayleigh- Jeans slope, becoming slightly steeper at high compactness, as discussed in 
SectionlOl 

The decomposition of the output spectrum into the components emitted at different radii is shown in Figure[l5] The contribution 
from large radius is in the optically thin regime, and has a harder spectrum than the (dominant) contribution near the photosphere. 
The overall normalization of the spectrum remains essentially constant, F Vm ~ Fq, at v em ~ uq when r es (r, , /i m i„) > 1 . 

6.5. Low-Frequency Spectral Slope 

Here we examine in more detail the effect that radiation transfer near a relativistic photosphere will have on the low-frequency 
spectrum. A Rayleigh- Jeans spectral slope arises from si de s cattering a monochromatic source in a locally spherical outflow, and 
is therefore maintained for a black body source (Section |63]l. The optical depth of low-frequency photons scales as T es (r, f? em ) oc 
oc r _2 i/7m and emission time t oc for photons emitted off the axis to the observer. The increase in optical depth also 
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then result from further upscattering by a process not considered in this paper. 
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FIG. 14. — Left panel: Slope d\n(L„)/d\ni/ of the spectra displayed in Figure [l2](blackbody source), showing a flat spectrum just below the scattered peak, 
steepening to a Rayleigh- Jeans slope at low frequencies. Solid lines: photon source attenuated by scattering added to the scattered spectrum. Dotted line shows 
blackbody, shifted arbitrarily in peak frequency for comparison. Right panel: Same outflow, but now a GRB-like photon source. Dotted line shows source, shifted 
in peak frequency. See text for discussion. 
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pushes the off-axis photosphere out to a radius r T (v em ) oc T" 2 ^^, where the rate of scatterings in a volume ~ r\ is proportional to 
(1 - /3)(Tp)n-y oc T~ 2 r T (v em )~ l oc v em . In the optically thin regime, one therefore finds a low-frequency spectrum F Vm oc t~ l oc v em , 
which hardens to F Vm oc v 2 m if every frequency is emitted from its photosphere. 

It is also worth examining briefly how this result for a steady outflow would be modified in the case where the outflow is 
impulsive. Consider a slightly simpler situation in which the seed photons flow radially (and are monochromatic). Then the 
frequency of the outgoing photon is purely a function of scattering angle, and we can consider the number of photons scattered 
the bulk frame of the magnetofluid in a time interval dt', 



d 2 N~, 



(70) 



Here n' ~ n^/T is the bulk frame photon density, and the differential scattering cross section varies mildly with the rest-frame 



scattering angle. Since v em /va = (l-/3)/(l-/3/i em ) = T 2 (l - (3)(l +/3/4 m ), we have 



" 7 o c 



d 2 N~, 



dv em dt> T 2 P{\-P)d^ 



I 

cm 



(71) 



The observed arrival time of a photon depends on emission angle and therefore frequency, dt^ = T(l - (3p em )dt' = v~mT(l - 
P)dt' . But integrating over the entire history of a pulse, the total number of photons emitted is flW 7 /aflni/ em oc v em . At very low 
frequencies, one must compensate for the expanded photosphere and reduced scattering rate, as outlined above. 



6.6. Semi-analytic Approximation to the Spectrum 

In the Monte Carlo evaluation of the scattered spectrum, we tested a simplified evaluation of the optical depth integral, using the 
small-angle approximation and assume linear growth of the Lorentz factor, T(r2) = (f2/r)T(f) ^> 1. Then equation (l65l l becomes 



7"csv : pern) - 



o T dN 
6cr dft 



_J_ , q2 
r 2 (r) em 



(72) 



The output spectra are hard to distinguish from those displayed in Figures ( [Tol l. ( fTTl i. 

It is also useful to work out the spectrum of singly-scattered photons in the case where the seed photons see a small to modest 
optical depth Since the matter is cold, it is simplest first to transform into its rest frame and consider the power scattered into 
solid angle dfl' em , 



d 3 E' 



dv'dSl'dt' 167r 



^ / c 



l+ik'-k'Y 



dtt'. 



(73) 



The rest-frame frequency and spectral intensity are v 1 = vq/T(\ + ftp') and l' v , = 7„/[T(l + (3p')] 3 . One sees in Figure |9]that the 
radiation field flows both forward and backward in the bulk frame close to the engine. At large distances the radiation field 
continues to collimate even as the V saturates, and V , is concentrated in the forward (anti-radial) direction. 
The power of the scattered radiation from a single charge is 



d 2 E 

dVemdt 



d 3 E' 



r dv> em dn> em dt> 
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dSl'dSl' 



1 + ^ 7 



(74) 



Using 



4 



l + (k'-k' em ) 2 =2n 2 3 + 3p' p' e 



-p 



I 2 
" Mem 



(75) 



gives 



d 2 E 3tt(Tt Iqvq _ 



dv em dt 8 f3T 4 



dp 



/ 3 + 3 / i Mem ZE Zl^B 
(1 + Ppl) 2 



(76) 



where p' em = p' em (p' , v em ) from equation ( |62l . The range of integration over pi is restricted if v em > 1. Since 1 + j3p' = (1 + 
/^MemV^em, we have (1 -j3)/v em < l+/3p' < (1 + /3)/j> em , and more generally 



1-/3 



1 



The flux of scattered radiation measured at a large distance is 

1 



< 1 + /V < (l+/3)min( 1, — 



dr dN d 2 E 
j3{r)c dft dv em dt 



(77) 



(78) 
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Normalizing to the incident radiation flux Fq = ttIqvq/x 2 , gives 




-3,-3/2/2 /2 / 2 
3 + 3^ Mem ZE 

r 4 /3 2 (i+/3/i') 2 



(79) 



Following equation (l72l . the prefactor can be written as 



3 or 



= -r 2 (r,)r es (r S) l). 



(80) 



7. SUMMARY AND CONCLUSIONS 



We have considered a very luminous and strongly magnetized outflow outside its scattering photosphere. The outflow is 
accelerated to a high Lorentz factor by a combination of the Lorentz force (which acts in a cold MHD flow) and the radiation 
scattering force (dominant in thermal fireballs). A range of radiation intensities is considered, extending from an almost cold 
flow to one in which the radiation and magnetic Poynting fluxes are comparable. The calculations described in this first paper 
assume that the poloidal magnetic field is strictly monopolar, which results in a near degeneracy between magnetic pressure 
gradient and curvature forces. Similar solutions will obtain for any part of an MHD outflow in which the magnetic flux surfaces 
are unfavorably curved and the Lorentz force remains weak. The opposing case, corresponding to a flared jet that breaks out of a 
confining medium, is examined in detail in paper II. 

The radiation force dominates the acceleration if the compactness \ > a at the photosphere. This inequality is easily satisfied 
if the outflow is optically thick near the engine: one has x( r r) ~ 6T 2 (r T )a in an outflow with comparable Poynting and photon 
energy fluxes. Radiative driving is especially efficient beyond the fast critical point, even though the magnetic field dominates the 
inertia of the outflow. The solutions we obtain for high radiation intensities can easily be rescaled to an outflow with relativistic 
bulk motion at a displaced photosphere: then the flow profile T(r) is linear in the inner parts of the outflow at both large and small 
scattering depths. 

We have considered the imprint of bulk Compton scattering on a photon seed with an exponential, high-frequency spectral 
cutoff. The spectrum is strongly modified when the radiation compactness is low enough that outward acceleration is dominated 
by the Lorentz force and the seed photon beam is wider than the Lorentz cone of the magnetofluid. In this situation, the mag- 
netofluid is pushed quickly to a high Lorentz factor outside its photosphere, where it feels a strong photon drag. Then the output 
spectrum extends above the seed peak frequency, with its low-frequency part depending on the shape of the seed. In the case of 
a blackbody seed, the spectral slope in between the seed thermal peak and the scattered peak is softer than Rayleigh- Jeans, but 
harder than is typical of GRBs at low frequencies. 

We have also considered a seed spectrum F v ~ const x e ~ hu l kT ° , representing the low-frequency part of the Band function absent 
the high-frequency tail. In that case, the scattered spectrum extends the flat portion of the seed spectrum upward in frequency. 
For both types of seed spectrum, the residual amplitude of the seed thermal peak that persists in the transmitted spectrum depends 
on the optical depth. A seed photon beam that is much wider than the Lorentz cone of the magnetofluid sees a large optical depth 
(compared with the optical depth of the more strongly beamed scattered photons), which means that the transmitted thermal peak 
is relatively weak. In the context of GRBs, the angular broadening of t he seed photons c ould result from scattering by a second, 
slower component of the outflow that is swept up at the head of the jet ([Thompson 2006)). 

In principle, no fine tuning of the optical depth surface is needed to make radiative and MHD acceleration competitive near 
the photosphere. The effective magnetization is much reduced below the photosphere, where the stress-energy of the photons 
couples to the matter and thence to the magnetic field: a e g ~ L P /L 7 . Once the radiation begins to stream freely, the outflow 
experiences both a rapid increase in magnetization, and a strong outward force from the self-collimating radiation field. This 
effect is examined in detail in Paper II. 



At first sight, one might associate the high-energy tail of the scattered photon spectrum with the observed high-energy tails of 
GRBs, and the seed photon energy with the observed spectral peak energy E p ^. But the calculated spectrum is relatively hard 
compared with the high-energy tails of GRBs, and it is limited in spectral width. 

Instead it appears more promising to identify the high-energy peak of the scattered spectrum with the measured E p ^ - at least 
in some bursts or possibly some phases of the burst emission. Then an additional source of dissipation, which is left out of our 
calculations, is needed to generate the high-energy tail. The seed thermal photons, generated deep in the outflow, provide a buffer 
that suppresses bursts with low E p ^, but upscattering allows a range of higher E p ^ values. We leave open here the nature and origin 
of the low-frequency seed spectrum, except to say that a hard Rayleigh- Jeans slope is by no means guaranteed if thermalization 
occurs at an intermediate optical depth. 

In this context, it is worth recalling some features of the Ama ti et ail (|2002) relation between £ p k an d the apparent is otropic 
energy £i S0 of GRBs. This relation is suggestive of jet breakout from the core of a Wolf-Rayet star (e.g. Thompson 2006), but it 
appears to represent a boundary in the E p ^-Ei so plane. The observed bursts have E p ^ lying on or above the Amati et al. line, with 
a strong deficit mainly below the line: for example, a signi ficant proportion of BATSE burst spectra are too hard to be consistent 
with this relation, independent of the (unknown) redshift (Nakar & Piran 2005). 

Evidence for a low frequency blackbody component has been found in some gamma-ray bursts (e.g. Axels son et al.l2012b . It is 
natural to try to accomodate such components in the present calculation, where the thermal seed should be partially transmitted. 
The measured soft blackbody component sits a factor ~ 1/30 below the burst peak, which requires a large broadening of the seed 
photons (by a factor ~ r# see d ~ 5 in angle, where 6> see d is the angular width of the seed component). If the decrease in the peak 
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frequency during the burst evolution were due mainly to a drop in r# see d, then one would expect to see a strong re-emergence 
of the seed thermal peak alongside the scattered peak, which is generally not observed. When considering the emission from 
relativistic outflows, one must always keep in mind the basic degeneracy between temporal and spectral degrees of freedom, and 
the angular degree of freedom. Therefore the presence of a soft thermal component in the spectrum may signal the presence of 
off-axis emission that does not, necessarily, interact with the emission zone of the harder component. 

Finally, onee should keep in mind that the scattered thermal radiation emitted by an electromagnetic outflo w could represent 
a subcomponent of the burst emission, with some other non-thermal (e.g. synchrotron) process dominating dPe'er et al.ll2006t 
Zhang & Yan 201 1). In that case, the photon field would have a smaller (but not necessarily negligible) influence on the outflow 
acceleration. 

7.2. Alternative Acceleration Mechanisms for Strictly Radial Magnetized Outflows 

We find that a hot electromagnetic outflow typically experiences a strong radiation force before it expa nds far enough tha t 
radial inhomogeneities become important. The acceleration of a static, bounded magnetic slab, studied by iGranot et al.l d201 ll) . 
could be relevant for the later stages of impulsive GRB outflows, outside a radius ~ cAf. But if a magnetized shell already 
moves relativistically at this radius, only a thin outer layer, comprising a fraction ~ 1 /2T 2 of the shell, would experience a strong 
outward magnetic pressure gradient force. Its interaction with slower material, swept up from a Wolf-Rayet star or a preceding 
neutron-rich wind, must then be taken into account. 

Other acceleration mechanisms have been suggested which depend on more complicated, non- ideal MHD effects, such a s 
the creation of a (net) outward pressure gradient force by reconnection of a toroidal magnetic field dDrenkhahn & Spruitli2 002). 
Zone s of alternating B^, separated by current sheets, are indeed present in the force-free solution to the oblique rotator ( Spitkovsky 
2006), in a zone straddling the rota tional equator. A n active dynamo operating in a GRB engine could also lead to stochastic 
reversals in the wind magnetic field (Thompson 2006). But when the increase in flow inertia associated with particle heating and 
radiation is taken into account, we have argued that cancelling even half the magnetic flux leads only to mildly relativistic radial 
motion. Magnetic reconnection plays a more natural role in the GRB phenomenon by modifying the gamma-ray spectrum, via 
particle heating and stochastic bulk motions. 

We thank the NSERC of Canada for financial support, and the referee for comments. 

APPENDIX 
ROTATING EMISSION SURFACE 

The photon source rotates rapidly in some cases, e.g. a rapidly rotating star such as a millisecond magnetar, or the merged 
remnant of a white dwarf binary. We can approximate the effect of a rotating emission surface by setting 




(Al) 



in equations d28l ), and d29i >. Here (3^^ is a constant representing the aberration of the outflowing photons at r = r s (x = 1). In 
this situation, plasma near the emission surface can more easily co-rotate with the radiation field while still being accelerated 
outward. 

The value of depends on the type of source. One has fl^R ~ Q.r s Jc = uj when the photons flow from the surface of a star 
of radius r s through a transparent wind. On the other hand, if the outflow is optically thick in a narrow radial zone close to the 
engine, then one expects (3^^ ~ (flr s /c)~ l ~ u~ l based on the conservation of angular momentum from the light cylinder out to 
the transparency surface (x = 1). 

In a first approximation, rotation of the photon source makes only small changes to the profiles of Lorentz factor and angular 
momentum in the outflow. For completeness, we discuss some of the detailed changes that do result in the flow parameters at 
the fast point. These effects are largest in intense radi ation fields with compactness \s > 4ct 4//3 /oj 2 . Then, as in our previous 
calculations with a non-rotating photon source (Section l4~4l i, the location of the critical point is determined by setting r c = r eq (r) 
[equation d32bl. Corotation of the fluid and the radiation field implies 




(A2) 



where x c ~ r c /3'/ 4 . When 0$ / 0, the fast speed ( |50l becomes 




(A3) 



At large compactness, the Lorentz factor and angular momentum at the critical point are given by 
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FIG. 16. — Effect of rotation of the photon source (compactness Xs) on the properties of the fast critical point: radius x c and flow variables T, flj, and C. Solid 
line: /3^ = l/u)=0.1 [equation (Al}]; dotted line: /3^j; = 0. Flow magnetization <r = 10 3 . 



One sees that the inertia added by the radial magnetic field can be significantly reduced when /3s ~ 1 /xuj, that is when f3^ R ~ u> 1 . 

A comparison of flow profiles with, and without, rotation of the photon source is made in Figure [16] There is little change 
in the critical point radius, but T c remains close to cr 1 / 3 at all values of the compactness, since we have taken /3$j{ ~ a; -1 . The 
magnetofluid rotates more rapidly near the critical point at large values of Xs- 
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